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1 Introduction 

It has been argued that differences observed between measurements of some ffavor 
observables and the corresponding Standard Model (SM) predictions may be due to 
beyond the SM (BSM) physics affecting the neutral 5-meson mixing processes [HIS]- 
Although the most recent analysis seem to indicate that there are not large BSM con- 
tribution to neutral 5- meson mixing |3], the future will bring new twists, and precise 
calculations of the non-perturbative inputs parametrizing the mixing in the SM and 
beyond are necessary for a thorough understanding of quark flavor physics. Lattice 
QCD can provide those non-perturbative theoretical inputs from first principles and 
with errors at the few percent level. 

The use of lattice QCD techniques can also shed light on the tension observed 
between the measured value of Br{B — )■ tv) and SM predictions by providing 
precise values of fs- Unitarity Triangle fits are very sensitive to fs and processes 
with potential to reveal NP effects depend on or /s^, so any improvement in the 
decay constant calculations is very important. 

In the next Sections I will summarize the status of calculations of both 5-meson 
decay constants and non-perturbative inputs relevant for the study oi — B mixing 
using lattice QCD techniques. I will only discuss realistic calculations with a complete 
error budget, simulations at several lattice spacings (and an extrapolation to the 
continuum), and including three dynamical quark flavors [Nj = 2 + 1 simulations). 



2 Decay constants 

Simulating heavy quarks on the lattice implies having to deal with discretization 
errors entering in powers of the mass in lattice units, amq. These corrections are 
not neghgible at typical lattice spacings a. In particular, the h quark can not be 
simulated at its physical mass on present lattices even with improved actions, since, 
typically, ami, > 1. There are two approaches that have been used to solve this 
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problem: simulating the heavy quarks with effective theories (such as heavy-quark ef- 
fective theories or non-relativistic QCD), and performing relativistic simulations with 
improved actions but with smaller masses than the physical bottom mass (but the 
same order or larger than the charm mass) and then extrapolating those results up to 
the physical rrib. The first approach is being used by the Alpha (heavy-quark effective 
theory, HQET), HPQCD (non-relativistic QCD, NRQCD), FNAL/MILC (Fermilab 
action), and RBC/UKQCD (non-peturbatively relativistic heavy-quark action) Col- 
laborations. One of the dominant systematic errors in these calculations is the one 
associated with the use of an effective theory. The collaborations that follow a rela- 
tivistic approach for b quarks and whose results I will mention here are ETMC (tm 
action) and HPQCD (HISQ action). 

There have been three lattice Nf = 2 + 1 calculations of fs and Jb^ in the last 
two years by the HPQCD [5l [6] (with a relativistic and a non-relativistic approach 
respectively) and the FNAL/MILC [7] collaborations, which have reduced the error 
to the 2.5% level. The smallest error in [B] for fs is achieved by determining the 
ratio fssl Ib using a non-relativistic description (NRQCD) of the h quarks together 
with the determination in [5], which employs relativistic h quarks. By using the 
ratio fBs/fs the dominant systematic errors associated with the effective NRQCD 
description partially cancel, so this determination is nearly free of uncertainties due 
to an effective description of the b quark. The average values of fs and fs^ from 
these three calculations are [8] 

fe = 190.6(4.7)MeV, /s, = 227.6(5.0)MeV , and ^ = 1.201(17) . (1) 

jBd 

The direct comparison of the results in Eq. ([1]) with experiment is problematic 
due to the need of the value of the CKM matrix element \Vub\ (whose inclusive and 
exclusive determinations disagree at the 3a level) and the ~ 2a disagreement of 
BaBar ^ and Belle's [10] measurements. Nevertheless, Belle new measurement seems 
to alleviate the tension between theory and experiment. 

In the next two years there will be new results for 5— meson decay constants with 
Nf = 2 + 1 and Nf = 2 + 1 + 1 from lattice collaborations using both an effective 
theory description (RBC/UKQCD [H]) and a relativistic description (ETMC0), as 
well updates with considerably smaller errors from FNAL/MILC ^13j . 

2.1 Neutral 5— meson mixing 

The current status of Nf = 2 + 1 lattice calculations of the non-perturbative quantities 
parametrizing the mass differences between the heavy and the light mass eigenstates 
in both the B'^ and 5° systems, as well as the SU(3) breaking ratio ^, is summarized 

1 Their Nf = 2 results [H] agree within la with Eq. 
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HPQCD 


FNAL/MILC 


RBC/UKQCD 




1.258(33) 


1.27(6) 


1.13(12) 


BbJBb^ 


1.05(7) 


1.06(11) 




HPQCD: 


Ibs \I Bbs = 


-- 266(6) (17) MeV, 


Bb, = 1.33(6) 


HPQCD: 


Ib^^/Bb^ = 


216(9)(13) MeV, 


Bb, = 1.26(11) 



Table 1: 5— meson mixing parameters. ^ is defined as the ratio of tlie parameters 
in tlie second and tliird rows. In tlie case wliere tliere are two errors, tlie first one is 
statistical and the second one systematic. 



in Tab. [H The HPQCD [H] and FNAL/MILC [I5] collaborations use the same hght 
quark formulation, but a different description for the b quarks. The exploratory 
study by the RBC/UKQCD collaboration uses heavy quarks in the static limit [T6]. 
The average of the results in Tab. [1] for ^ gives the value ^ = 1.251 ± 0.032. The 
calculations whose results are shown in Tab. [T]were not optimized to extract the bag 
parameters but rather the matrix elements, so the errors for the bag parameters can 
be significantly reduced to about a 3% error in on-going calculations. 

There is not yet a finalized calculation of the matrix elements needed for the 
determination of the decay width differences, ATd^s, in the continuum limit and with 
Nf = 2 + 1 flavors of sea quarks, but preliminary results for the relevant matrix 
elements by FNAL/MILC can be found in [T7] . 

Beyond the SM the mixing parameters can have contributions from AB = 2 four- 
fermion operators which do not contribute in the SM. The matrix elements of the 
five operators in the complete basis describing AB = 2 processes, together with the 
Wilson coefficients for those operators calculated in a particular BSM theory and the 
experimental measurements of the mixing parameters, can provide very useful con- 
straints on that BSM theory. Again, there is not a final unquenched lattice calculation 
of the matrix elements of all the operators in the AB = 2 effective hamiltonian, but 
FNAL/MILC presented preliminary results for the complete basis in [17]. 

The authors of Ref. [18] suggested that the branching fractions of the rare decays 
Bg — )■ fi^fi~ (for q = s,d) could be determined from the experimental measurement of 
the mass difference in the neutral i?q-meson system, AMg, and the lattice calculation 
of the bag parameter Bb^ using 

AMg ^^'''^''''riBWMwsin^e^) ^ S{x,) Bg ^ ' 

In order to compare experimental measurements and the theory predictions for the 
decay rate of B^, one must include the effects of a non- vanishing AP^ [1^. This 
can be done in the SM by rescaling the theory prediction by 1/(1 — t/s), where ys = 
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rB,Ar,/2 [19]. Multiplying Eq. (ED by this factor for the 5° fi+fx' decay and 
using the HPQCD determination of the bag parameters Bb^ = 1.33(6) and Bb^ = 
1.26(11) [H]; together with tb, = 1.497(15)ps, r^, = 1.519(7)ps [20], and AF, = 
0.116(19)ps-^ [21], one gets 

Br{Bs Ai+/u")|y, = (3.65 ± 0.20) x 10"^ , 

Br{Bd^ fi^fi') = (1.04 ± 0.09) X 10^^°. (3) 

The direct calculation of these branching fractions has become competitive with 
the one in ([3]) [22j thanks to the recent improvements in the calculation of the B- 
meson decay constants on the lattice summarized in Sec. [2l Including the correction 
factor 1/(1 — Us) for the Bs — ?■ fi'^fi~ decay rate, and using the same inputs as in 
Ref. [22i except for /b and /b^, for which I use the averages described in Sec. [H and 
tb^, which I take equal to its PDG 2012 value, tb^ = 1.497(15)ps, I get 

Br{B, ^ V~) \ys = (3-64 ± 0.23) ■ 10"^ , 

Br{Bd^ fx+fi') = (1.07 ± 0.10) ■ 10"^°. (4) 

The agreement between the two set of numbers in and (jlj) is excellent. This gives 
us confidence in the SM prediction for these branching fractions, and this confidence 
will increase when we have results for the bag parameters entering in ([3]) from the 
on-going lattice calculations described above. This is very important since the LHC 
bounds are now approaching the SM predictions 0, especially for B^ decays, Br{Bs — t- 
I^^I^~)\lhc < 4.2 ■ 10-9 at 95% CL [21]. 

3 Outlook and future prospects 

In the near future there will be results for the S-meson mixing parameters and the 
decay constants at the percent level and with different lattice formulations for both 
light and heavy quarks. Some of the results expected for next year are: the first 
unquenched determinations in the continuum limit of matrix elements needed for the 
calculation of the decay width differences AF^ the analysis of B mixing in BSM 
theories, and the study of short-distance contributions to — D mixing. 

Further improvement will be achieved by using simulations at the physical u and 
d masses (some preliminary results for Jb are already available from the HPQCD 
collaboration), including the effect of dynamical charm quarks, and extending the 
use of relativistic actions to the analysis of B^ — B mixing. 

^ After this conference, the LHCb Collaboration announced the rst evidence for one of these two 
processes [23], Br{Bs — ^ /i^/i^) = (3.2^i['2) • lO"^. This result agrees with the SM prediction given 
in this proceeding. 
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